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Abstract 

Fas and its ligand (FasL) play an important role in apoptosis and carcinogenesis. Therefore, the potential association of 
polymorphisms in the Fas (-670A>G, rs1800682; -1377G>A, rs2234767) and FasL (-8440T, rs7631 10) with cancer risk has 
been widely investigated. However, all the currently available results are not always consistent. In this work, we performed a 
meta-analysis to further determine whether carriers of the polymorphisms in Fas and FasL of interest could confer an altered 
susceptibility to cancer. All relevant data were retrieved by PubMed and Web of Science, and 52 eligible studies were 
chosen for this meta-analysis. There was no association of the Fas -670A>G polymorphism with cancer risk in the pooled 
data. For the Fas -1377G>A and FasL -844C>T polymorphisms, results revealed that the homozygotes of -1377A and -844C 
were associated with elevated risk of cancer as a whole. Further stratified analysis indicated markedly increased risk for 
developing breast cancer, gastric cancer, and esophageal cancer, in particular in Asian population. We conclude that carriers 
of the Fos-1377A and the FasL -844C are more susceptible to the majority of cancers than non-carriers. 
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Introduction 

With new cases and mortality increased dramatically, cancer 
has become the major public health burden worldwide. For this 
reason, novel diagnostic markers are needed urgently for early 
detection and prevention of cancer. However, carcinogenesis is a 
complicated biological process that is not fully understood. It is 
generally believed that interactions of low-penetrance susceptibil- 
ity genes with environmental factors might contribute to carcino- 
genesis [1]. As one of the important low-penetrance genes, Fas is 
considered to be a potential cancer susceptibility gene. This is 
because Fas (TNFSF6, CD95, or APO-1) is a cell surface receptor 
involved in apoptotic signal transmission in many cell types and 
interacts with its natural ligand Fas ligand (also known as FasL) to 
initiate the death signal cascade that leads to apoptotic cell death 
[2,3]. Furthermore, in these two genes, there are several 
functionally significant polymorphisms, such as the — 670A>G 
and — 1377G>A in the Fas promoter region, and the — 8440T 
in the FasL promoter region, because they might be associated 
with cancer risk, including cervical cancer [4—9], gastric cancer 
[10-15], breast cancer [16-21], lung cancer [22-25] and so on. 
However, all available results are not always consistent with one 
another, partially because of the small sample size of some 
published studies, different ethnic backgrounds, publication bias, 
and little effect of the polymorphisms on cancer risk. Therefore, 
it's necessary to retrieve and pool all eligible data to further 



determine whether these genetic polymorphisms could be at 
increased risk for developing cancer and to what extent 
heterogeneity existed across all the studies. 

Materials and Methods 

Identification and eligibility of relevant studies 

Two online medical databases, PubMed, and Web of Science, 
were searched (updated February 2013), using the search terms 
"Fas/CD95/TNFSF6/APO-l", "FasL/CD95L", "polymor- 
phism/genetic variation" and "cancer/carcinoma/tumor"). The 
literature search was limited to English articles. In addition, more 
studies were also identified by manual search based on the 
references provided in the retrieved studies. The inclusion criteria 
were prespecified as below: (1) be a case-control study, (2) evaluate 
association between the Fas and/or FasL polymorphisms and 
cancer risk, (3) present sufficient data to calculate an odds ratio 
(OR) with 95% confidence interval (CI), and (4) list genotype 
frequency. Moreover, the studies without raw data, or those that 
were case-only studies, case reports, editorials, and review articles 
(including meta-analyses) were eliminated. 

Data extraction 

Information was extracted carefully from all eligible articles 
independently by two authors (Yeqiong Xu and Bangshun He) 
according to the above inclusion and exclusion criteria. Discrep- 
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1 68 of records identified through PubMed searching 



120 of records identified through Web of Knowledge 



288 Potentially relevant trials identified 







80 duplicates removed 


► 

r 


208 of re< 


;ords screened 





148 of full-text articles excluded (59 not 
the variants of interest , 6 other SNP, 
45 not cancer studies, 21 not 
case-control studies, 4 reviews, 6 meta, 
7 not English) 



1 

60 of full-text articls 
article 


r 

■s assessed for eligibility 
s assessed 






► 


8 articles without independent data 





52 of studies included in qualitative synthesis 



Figure 1. Flow chart of studies identified according to inclusion and exclusion criteria. 

doi:1 0.1 371 /journal. pone.0090090.g001 



ancies were resolved by extensive discussion in our research team. 
The characteristics of enrolled studies were extracted as below: the 
first author's last name, year of publication, country of subjects, 
ethnicity, type of cancer, the source of controls, genotyping 
method (whether PCR was performed using a dual-labelled 
TaqMan probe with a specific 3'base to detect the SNPs or 
whether an RFLP method was used), the number of matched cases 
and controls, polymorphism sites, and P value for Hardy- 
Weinberg equilibrium (HWE) as summarized in Table 1. 

Genotype-gene expression correlation analysis 

The International HapMap Project (http://hapmap.ncbi.nlm. 
nih.gov/) was used to obtain data of the Fas and FasL genotypes 
determined in 270 enrolled subjects. Meanwhile, the mRNA 
expression data of these enrolled subjects were available online 
from SNPexp (http://app3.titan.uio.no/biotools/help. 
php?app = snpexp) as described in the previous studies [26,27]. 
In brief, these data were obtained from the HapMap phase II 
release 23 data set consisting of 3.96 million SNP genotypes from 
270 subjects of three populations, including 90 European (CEU), 
90 Asian (45 Chinese, 45 Japanese), and 90 Yoruba (YRI) subjects 
[28] . Additionally, the mRNA expression data were derived from 
the lymphoblastic cell lines from the same 270 subjects [29]. 

Statistical analysis 

Crude ORs with 95% CIs were used to assess the strength of 
association between the polymorphisms in fiw-670A>G, Fas - 
1377G>A, and FasL -844T/C and cancer risk. The pooled ORs 
were estimated for dominant model (variant homozygotes + 
heterozygous vs homozygous reference), recessive model (variant 
homozygotes vs heterozygous + homozygous reference), homozy- 
gote comparison (variant homozygotes vs homozygous reference), 



heterozygote comparison (heterozygous vs homozygous reference) 
and allelic comparison in the polymorphisms, respectively. 
Stratified analyses were performed by the type of cancer (that 
with only one study was grouped together as 'other cancers'), 
ethnicity, source of controls and genotyping method. Heteroge- 
neity across the studies was evaluated by using the Chi-square test 
based Q_-statistic test, and it was considered statistically significant 
when Phetemgendty (P/,)<0-05. The data were combined using 
random-effects model (the DerSimonian and Laird method) [30] 
in the presence of heterogeneity (P<0.05 or / >50%), or fixed- 
effects model (the Mantel-Haenszel method) models [31] was 
chosen to use in the absence of heterogeneity (P>0.05 or 
I <50%). Moreover, sensitivity analyses were performed to assess 
the stability of the results. Publication bias was evaluated 
graphically by using funnel plots and statistically by the Egger's 
linear regression test. HWE of the three polymorphisms was 
assessed using a web-based program (http://ihg.gsf.de/cgi-bin/ 
hw/hwal.pl). All statistical tests were performed with ST ATA 
1 1.0 and SPSS 20.0. All the P values were two-sided. 

Results 

A total of 52 studies were enrolled in this meta-analysis 
(Figure 1). The major characteristics of the 52 selected studies are 
summarized in Table 1. The study carried out by Bye et al [32] 
analyzed individuals of African or Mixed ethnicity, and thus was 
divided into two studies. Similarly, the studies reported by Ho et al 
[33] and Ueda et al [34] investigated two and three types of 
cancer, and therefore, these two studies were cited as two studies 
and three studies, respectively (Table 1). 

For the Fas -670A>G polymorphism, there was no association 
in the pooled analysis. In the subgroup analysis, statistically 
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Table 2. Stratified analyses of the Fas -670A>G (rsl 800682) polymorphism and cancer. 





Variables 


n a 


GG+AG vs AA 






GG vs AG+AA 






G vs A 






OR(95%CI) 




I 2 


OR(95%CI) 




I 2 


OR(95%CI) 




I 2 


Total 


44 


1.01(0.94, 1.09) c 


<0.0001 


47 '.1 


1.04(0.96, 1.1 2) c 


0.003 


40.9 


1. 02(0.97,1. 06) c 


0.005 


39.4 


Cancer type 


Cervical cancer 


9 


1. 05(0.79,1. 40) c 


-C0.0001 


74.5 


0.92(0.69, 1.22) c 


0.006 


62.8 


0.99(0.86,1. 14) c 


0.013 


58.5 


Gastric cancer 


5 


1.08(0.91,1.28) 


0.340 


11.6 


0.97(0.79,1.21) 


0.978 


0.0 


1.03(0.91,1.15) 


0.735 


0.0 


Esophageal cancer 


4 


1.02(0.85,1.21) 


0.459 


0.0 


1.21 (0.86, 1.69) c 


0.017 


70.4 


1.10(0.99,1.23) 


0.215 


32.9 


Breast cancer 


4 


1.01(0.90,1.14) 


0.325 


13.4 


1.03(0.90,1.18) 


0.062 


59.1 


1.02(0.94,1.10) 


0.259 


25.5 


Prostate cancer 


3 


0.83(0.70,0.98) 


0.155 


46.4 


0.82(0.66,1.01) 


0.346 


5.8 


0.87(0.77,0.97) 


0.163 


44.8 


Ovarian cancer 


2 


0.87(0.66,1.15) 


0.952 


0.0 


0.85(0.57,1.28) 


0.622 


0.0 


0.90(0.74,1.09) 


0.745 


0.0 


Bladder cancer 


2 


1.01(0.77,1.33) 


0.588 


0.0 


1. 00(0.47,2. 16) c 


0.043 


75.6 


1.03(0.85,1.24) 


0.491 


0.0 


Skin cancer 


2 


1.08(0.91,1.27) 


0.414 


0.0 


1.02(0.86,1.23) 


0.483 


0.0 


1.04(0.93,1.16) 


0.902 


0.0 


Nasopharyneal cancer 


2 


1.55(0.75,3.24) c 


0.017 


82.4 


1.39(0.69,2.79) c 


0.042 


75.8 


1. 33(0.80,2.1 9) c 


0.008 


85.6 


Melanoma 


2 


0.79(0.64,0.97) 


0.765 


0.0 


0.96(0.77,1.21) 


0.790 


0.0 


0.90(0.78,1.02) 


0.725 


0.0 


Lung cancer 


2 


0.82(0.65,1.04) 


0.852 


0.0 


1.07(0.82,1.40) 


0.906 


0.0 


0.94(0.81,1.10) 


0.984 


0.0 


Other cancers 


7 


1.08(0.96,1.22) 


0.373 


7.3 


1.15(0.99,1.32) 


0.747 


0.0 


1.08(1.00,1.17) 


0.528 


0.0 


Ethnicity 


Asian 


25 


0.97(0.88,1. 06) c 


0.004 


48.3 


1.01(0.89, 1.1 5) c 


0.003 


49.3 


0.99(0.93,1. 05) c 


0.030 


37.8 


Caucasian 


13 


1.03(0.95, 1.12) 


0.120 


32.8 


1.00(0.92, 1.09) 


0.277 


16.5 


1.01(0.96,1.06) 


0.277 


16.6 


African 


3 


1.72(1.24,2.38) 


0.288 


19.6 


1.23(0.78,1.95) c 


0.039 


69.1 


1. 25(0.90,1. 74) c 


0.022 


73.9 


Mixed 


3 


1.10(0.82, 1.48) 


0.607 


0.0 


1.28(0.99, 1.65) 


0.803 


0.0 


1.15(0.97,1.37) 


0.610 


0.0 



a Number of comparisons. 

h P value of Q-test for heterogeneity test. 

c Random-effect model was applied when P value for heterogeneity < 0.05; otherwise, fixed-effect model was applied. 
Statistically significant results were in bold. 
doi:1 0.1 371 /journal.pone.0090090.t002 



significantly decreased risk was observed in prostate cancer and 
melanoma for GG+AG vs AA comparison model, whereas there 
was significandy increased risk among those of African ancestry for 
GG+AG vs AA models (all data shown in Table 2). 

For Fas -1377G>A polymorphism, significantly increased 
cancer risks were observed in AA vs GG (Figure 2) and AA vs 
GA+GG comparison models in the overall analysis. In the 
subgroup analysis by cancer type, a significantly increased risk 
was observed in breast cancer for all comparison models. 
Meanwhile, increased risks were found for the comparison of 
AA vs GG and AA vs GA+GG in gastric cancer and esophageal 
cancer. In addition, a borderline decreased cancer risk was found 
in melanoma for GA vs GG and AA+GA vs GG comparison 
models (all data shown in Table 3). 

For FasL -8440T polymorphism, significandy increased 
cancer risks were observed in CC vs TT (Figure 3), CC+CT vs 
TT and CC vs CT+TT in the overall analysis. When the analysis 
was stratified by genotyping method, an increased cancer risk was 
observed in studies carried out by PCR-RFLP (shown in Table 4). 

Overall effects for alleles 

Allele comparisons were also conducted in the meta-analysis. 
However, no significant associations were found in Fas -670A>G 
polymorphism and cancer risks (shown in Table 2). 

There was borderline association between Fas -1377G>A 
polymorphism and cancer risks for A allele vs G allele in the 
overall analysis. In the subgroup analysis by cancer type, opposite 
results were shown between breast cancer and melanoma (shown 
in Table 3). 



For FasL -8440T polymorphism, in the subgroup analysis of 
genotyping method, an increased cancer risk was found in die 
studies carried out by PCR-RFLP (shown in Table 4). 

The Fas and FasL mRNA expression by genotypes and 
population 

The Fas and FasL mRNA expression levels were stratified by 
genotype (shown in Table 5) and population (shown in Table 6) 
groups. In the genotype subgroup analysis, significant associations 
between mRNA expression levels and Fas -670A>G were 
observed in all populations (GA: P= 0.043), especially in Asian 
population (GG: P= 0.0003; dominant: P= 0.003; recessive: 
P= 0.001). Meanwhile, significant differences between mRNA 
expression levels and FasL -8440T were observed in Asian 
population (recessive: P= 0.001). In the population-subgroup 
analysis, decreased expression of Fas was found in YRI (Yoruba 
in Ibadan) population than in the CEU population (P= 0.002). 

Test of heterogeneity 

There was significant heterogeneity across the studies focused 
on these three polymorphisms as evaluated by Q_-test. Then, we 
evaluated the heterogeneity for dominant model comparison by 
subgroups (cancer type, ethnicity, source of controls and genotyp- 
ing method). As a result, ethnicity {x 2 = 13.44, degree of freedom 
= 3, P h = 0.004) and cancer type (x Z = 22.26, degree of freedom 
= 11, P h = 0.022), but not source of controls [yf= 1.49, degree of 
freedom = 1, P h = 0.222) or genotyping method (x = 1 .48, degree 
of freedom =4, P h = 0.830) contributed to substantial heteroge- 
neity of the Fas -670A>G polymorphism. For the Fas -1377G>A 
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Figure 2. Forest plots of effect estimates for fas -1377G>A polymorphism (AA vs GG). For each of the studies, the estimation of OR and its 
95% CI is plotted with a box and a horizontal line. Filled diamond pooled OR and its 95% CI. 
doi:1 0.1 371 /journal.pone.0090090.g002 



polymorphism, the test revealed cancer type (% = 22.60, degree of 
freedom = 8, P h = 0.004), but not ethnicity (^ = 4.81, degree of 
freedom = 3, = 0.187), source of controls (x 2 = 0A2, degree of 
freedom = 1, P h = 0.518), or genotyping method (x = 0.51, degree 



of freedom = 3, P/,— 0.917) contributed to substantial heterogene- 
ity. For the FasL -844C>T polymorphism, genotyping method 
(X =9.21, degree of freedom = 3, P/, = 0.027), but not cancer type 
(X 2 = 4.33, degree of freedom = 7, P h = 0.741), ethnicity (/ = 5.64, 
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Figure 3. Forest plots of effect estimates for fa5i-844C>T polymorphism (CC vs TT). For each of the studies, the estimate of OR and its 
95% CI is plotted with a box and a horizontal line. Filled diamond pooled OR and its 95% CI. 
doi:1 0.1 371 /journal.pone.0090090.g003 



degree of freedom = 3, P/,= 0.131), or source of controls 
(X 2 = 0.08, degree of freedom =1, = 0.777) contributed to 
substantial heterogeneity. 

Sensitivity analyses 

To assess the stability of the results and the source of the 
heterogeneity, sensitivity analysis was performed by sequential 



removal of each individual eligible study. For Fas -670A>G and 
FasL -8440T polymorphisms, statistically similar results were 
observed after sequential removal of individual study in dominant 
and homozygote model, respectively, and the summary ORs in the 
other genetic models were not materially altered, suggesting that 
the results were stable. For the Fas -1377G>A polymorphism, 
sensitivity analysis indicated that study by Shao et al [38] was 
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responsible for heterogeneity. The heterogeneity was decreased 
when this study was removed (AA+GA vs GG: P h = 0.075, 
/ = 26.5). Although the genotype distribution in 1 1 studies (listed 
in Table 1) didn't follow HWE, the corresponding summary ORs 
were not materially altered with or without including these studies 
for the three polymorphisms. In addition, no other single study 
altered the pooled ORs by sensitivity analysis. 

Publication bias 

To assess the publication bias, Begg's funnel plot and Egger's 
test were performed and the shapes of funnel plots didn't show any 
obvious asymmetry in all genetic models of the three polymor- 
phisms (Figure 4A-C). Therefore, to provide statistical evidence of 
funnel plot symmetry, Egger's test was performed for each of these 
polymorphisms and the results confirmed the absence of 
publication bias (P>0.05). 

Discussion 

Fas, a potent member of the death receptor family, plays a 
crucial role in apoptotic signaling in many cell types [40]. 
Meanwhile, interactions between Fas and its receptor FasL trigger 
the death signal cascade, and subsequently induce apoptotic cell 
death [41]. Previous studies have indicated that down-regulation 
of Fas expression and/ or up-regulation of FasL expression could 
be detected in many types of human tumors [42,43] . The reason 
may be that down-regulation of Fas could protect tumor cells from 
elimination by anti-tumor immune responses, whereas up-regula- 
tion of FasL could increase the ability of tumor cells to 
counterattack the immune system by inducing apoptosis 
~v [44,45,46]. Therefore, it is believed that Fas and FasL play a 

o. crucial role in carcinogenesis. Given the important roles of Fas and 

% FasL in carcinogenesis process, it is biologically plausible that Fas 

and FasL polymorphisms that possess the potential to influence the 

0 expression of Fas and/ or FasL may be associated with cancer risk. 
„ Therefore, associations between the Fas -670A>G, Fas - 

u 

it 1377G>A and FasL -8440T polymorphisms and cancer risk 

■6 were determined in this meta-analysis. 

£ In this meta-analysis, 52 published studies were enrolled to 

.h determine the association between the three potentially functional 

1 polymorphisms within the Fas and FasL and cancer risk. This study 
I revealed that the Fas -1377G>A and FasL -8440T, but not the 
q Fas -670A>G polymorphisms were associated with significandy 
V increased overall cancer risk. Previous studies have identified that 
.|T the -1377A allele had markedly reduced ability to bind transcrip- 
§ tion factor stimulatory protein 1 as compared with the -1377G 
£ allele, whereas the -670A and G alleles had similar ability to bind 
|i transcription factor signal transducers and activators of transcrip- 
£ tion 1 (STAT 1 ) [47] . As the Fas - 1 3 7 7 A allele reduced the ability to 
^ bind transcription factor stimulatory protein 1 that is a crucial 
> transcriptional activator, the expression of Fas was decreased in 
c ^ carriers of the Fas -1377AA genotype as expected, but the Fas - 
■f g 670G allele didn't influence the expression of Fas [47,48]. 
"a o ° Therefore, it is reasonable that the Fas -1377A allele increased 
S- g o the overall cancer risk, and that the Fas -670G allele had no 

marked effect on overall cancer risk, which was consistent with our 
results. For the FasL -844T>C polymorphism, which is located in 
a binding motif for transcription factor CAAT/ enhancer binding 



O ™ 41! CTi 



° ^ <D P O 
O O 5 

S g « E | S protein P, could influence the promoter activity of the FasL gene 
£ ° 6,1 ? [49] . Additionally, it has been proposed that compared with the - 

£ J ° | = s 844T allele, -844C allele strongly increased the expression of FasL 

~ -i -o « 2 on T cells and was associated with an enhanced rate of activation- 



E 

~ ii 



I S '6 induced cell death of T cells, which may lead to less powerful 



immune surveillance and increase the susceptibility to cancer [6] . 
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Table 5. Fas and FasL mRNA expression by the genotypes of SNPs, using data from the HapMap 1 . 





Fas -670A>G 










FasL -8440T 








Population 


Genotypes 


No. 


Mean ± SD 


F* 


Ethnicity Genotypes 


No. 


Mean ± SD 


P* 


CEU 3 


AA 


23 


8.79 ±0.36 




CEU 3 CC 


76 


5.94±0.07 






GA 


46 


8.87 ±0.28 


0.321 


CT 


5 


5.89±0.07 


0.137 




GG 


12 


8.74±0.36 


0.687 


TT 


0 


— 


— 




Dominant 


58 


8.84±0.30 


0.511 


Dominant 


5 


5.89±0.07 


0.137 




Recessive 


69 


8.84±0.31 


0.292 


Recessive 


81 








YRI 3 


AA 


6 


8.58±0.33 




YRI 3 CC 


0 


— 


— 




GA 


25 


8.70 ±0.31 


0.402 


CT 


28 


5.94±0.06 


— 




GG 


53 


8.67 ±0.30 


0.450 


TT 


53 


5.95 ±0.06 


— 




Dominant 


78 


8.58±0.33 


0.410 


Dominant 


81 


5.95 ±0.06 


— 




Recessive 


31 


8.67 ±0.31 


0.987 


Recessive 


28 


5.94±0.06 


0.493 


Asian 3 


AA 


28 


8.65 ±0.29 




Asian 3 CC 


0 





— 




GA 


36 


8.78 ±0.26 


0.059 


CT 


50 


5.96±0.06 


— 




GG 


21 


8.98 ±0.30 


0.0003 


TT 


33 


5.91 ±0.06 






Dominant 


57 


8.85 ±0.29 


0.003 


Dominant 


83 


5.94±0.06 






Recessive 


64 


8.72 ±0.28 


0.001 


Recessive 


50 


5.96±0.06 


0.001 


All 3 


AA 


57 


8.70±0.33 




All 3 CC 


76 


5.94±0.07 






GA 


107 


8.80±0.28 


0.043 


CT 


83 


5.95 ±0.06 


0.163 




GG 


86 


8.76±0.33 


0.297 


TT 


86 


5.94±0.06 


0.913 




Dominant 


193 


8.78±0.30 


0.081 


Dominant 


169 


5.95 ±0.06 


0.390 




Recessive 


164 


8.76 ±0.30 


0.871 


Recessive 


159 


5.95 ±0.06 


-C0.0001 



CEU: 90 Utah residents with ancestry from northern and western Europe; YRI: 90 Yoruba in Ibadan, Nigeria; Asian: 45 unrelated Han Chinese in Beijing and 45 unrelated 
Japanese in Tokyo. 

'Genotyping data and mRNA expression levels for Fas and FasL by genotypes were obtained from the HapMap phase II release 23 data from EBV-transformed 

lymphoblastoid cell lines from 270 individuals. 

2 Two-side Student's f test within the stratum was used. 

3 There were missing data for unavailable genotyping data. 

Statistically significant results were in bold. 

doi:1 0.1 371 /joumal.pone.0090090.t005 



The Fas -670GG genotype was associated with decreased risk of 
prostate cancer and melanoma according to the cancer type 
subgroup analysis. It was suggested that Fas -670A>G polymor- 
phism might have the same effect on these two cancers. However, 
these results were based on 44 studies, which could affect the 
results owing to small amount of studies. Therefore, to draw a 
more precise conclusion, more related studies are needed. 

For the Fas -1377G>A polymorphism, this study revealed that 
those who carried the -1377AA genotype had an increased risk for 
breast cancer, gastric cancer and esophageal cancer, while the 
melanoma risk was decreased. As described above, the different 
risk factors could contribute to the discrepancies. Also other 
unidentified causal genes would influence the effect of this 
polymorphism on different cancers. 

For the FasL -8440T polymorphism, the -844CC associated 
with increased cancer risk was observed in gastric cancer, 
esophageal cancer, and ovarian cancer among the previous 
studies, indicating that this polymorphism had similar effect on 
these three cancers. Although these cancers had different 
mechanisms of carcinogenesis, small amount of studies, 
publication bias, and other unidentified causal genes would be 
the result of the discrepancies, which contributed to the similar 
association between the FasL -8440T polymorphism and three 
cancers. 

In the subgroup analysis by ethnicity, an increased cancer risk in 
carriers of the Fas -670GG genotype was found in African, while 



the result of mRNA expression showed that GG genotype 
expressed higher levels of Fas in Asian populations. Meanwhile, 
the previous studies showed increased cancer risk in carriers of the 
Fas -1377AA and FasL -844CC genotype were found in Asian 
subjects, which was evidenced in mRNA expression by genotypes 

Table 6. Fas and FasL mRNA expression by the ethnicity, 
using data from the HapMap 1 . 



Fas -670A>G FasL -8440T 

Ethnicity No. Mean ± SDP 2 Ethnicity No. Mean ± SDP 2 

CEU 3 81 8.83±0.31 CEU 3 81 5.94±0.07 

YRI 3 84 8.67±0.30 0.002 YRI 3 81 5.95 ±0.06 0.120 

Asian 3 85 8.79±0.30 0.391 Asian 3 83 5.94±0.06 0.398 

CEU: 90 Utah residents with ancestry from northern and western Europe; YRI: 90 
Yoruba in Ibadan, Nigeria; Asian: 45 unrelated Han Chinese in Beijing and 45 
unrelated Japanese in Tokyo. 

Genotyping data and mRNA expression levels for Fas and FasL by genotypes 

were obtained from the HapMap phase II release 23 data from EBV-transformed 

lymphoblastoid cell lines from 270 individuals. 

2 Two-side Student's f test within the stratum was used. 

3 There were missing data for unavailable genotyping data. 

Statistically significant results were in bold. 

doi:1 0.1 371 /joumal.pone.0090090.t006 
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Figure 4. Begg's funnel plot of Egger's test for publication bias for three polymorphisms. Each circle represents as an independent study 
for the indicated association. Log[OR], natural logarithm of OR. Horizontal lines mean effect size. A: Begg's funnel plot of publication bias test for Fas 
-670A>G polymorphism. B: Begg's funnel plot of publication bias test for Fas -1377G>A polymorphism. C: Begg's funnel plot of publication bias 
test for FasL -8440T polymorphism. 
doi:1 0.1 371 /journal.pone.0090090.g004 



in Asian populations. However, this association was not proved in 
other ethnicities. The discrepancies in racial backgrounds and 
environment they lived in would lead to the differences. In 
addition, these polymorphisms might be masked by the presence 
of other unidentified causal genes involved in carcinogenesis. Due 
to the small size of population for the ethnicities, well-designed, 
large randomized case-control studies should be performed. 

The pooled results of this study may be affected by polymor- 
phism genotyping methods applied in the enrolled studies. 
Previous studies revealed that the pooled results of the Fas - 
670A>G polymorphism were not affected by the studies with 
genotyping methods of both PCR-RFLP and TaqMan. While Fas 
-1377AA genotype carriers increased cancer risk in the studies 
using PCR-RFLP but not TaqMan, and similar result was found 
in the FasL -844CC genotype carrier. The discrepancy across the 
studies applied different polymorphism genotyping methods may 
result from the different sensitivity and accuracy of genotyping 
methods. Meanwhile, the quality control is crucial to cause 
discrepancy as well. In general, studies [12,17,50] selected 10% 
repeated, random sample of subjects to test twice by standard 
genotyping method or different investigator, which was used to 
confirm the accuracy of results, while Mandal et al [51] and Ter- 
Minassian et al [22] tested 5% repeated samples. As a result, the 
consistency rate of quality control was 100% in almost all studies. 
However, the study by Crew et al [19] showed that the consistency 
rate was 100% for Fas -1377G>A, 94% for Fas -670OA and 
96% for FasL -8440T. Therefore, the results of further studies 
should be confirmed by a standardized genotyping method. In 
addition, the limited amount of studies would also contribute to 
the discrepancy. 

Heterogeneity is an important factor which can interpret the 
results of the meta-analysis. Therefore, we stratified the studies by 
cancer type, ethnicity, source of controls and genotyping method, 
respectively. The results showed that the main heterogeneity 
existed for cancer type and ethnicity. The reason might be that 
different cancers have different mechanisms of carcinogenesis. 
Virus infections, hormone levels, smoking, drinking, family history 
all could contribute to the different cancers. Meanwhile, different 
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